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Dispositif et precede de d ecodage rnhuste de codes arithmetiq ues 

L'invention concerne la compression/decompression de donnees numeriques, notamment 
5 pour des signaux multimedia (audio, image, video, parole), et la transmission robuste de ces 
donnees sur des reseaux bruites, tels que des reseaux sans fils, des communications mobiles. 

Pour reduire le debit de transmission des donn6es numeriques, on les comprime, en 
cherchant a s'approcher du maximum theorique que les specialistes appellent 'Tentropie du 
1 0 signal". Pour cela, on utilise souvent des codes statistiques aussi appeles codes a longueurs 
variables, par exemple les codes de Huf&nan. Ces codes ontneanmoinsl'inconv6nientd'etre 
tres sensibles aux erreurs de transmission. L'inversion de un bit peut conduire a une 
desynchronisation du d6codeur, qui a pour consequence un decodage errong de toutes les 
donn6es qui suivent la position du bit errone. 

15 

Les solutions actuelles de compression, transmission et decompression de signaui 
multimedia sur reseau, sur lesquelles on reviendra, se fondent sur l'hypothese qu'urie 
certaine qualite de service du transport des donnees est garantie. En d'autres termes, ell6s 
supposent que les couches de transport et de liaison, en s'appuyant sur l'utflisation de codes 

2 0 correcteurs, vontpermettre d'atteindre un taux d'erreur residuel quasi nul (c'est-a-dire vu par 
1' application de compression et de decompression). Mais cette hypothese de taux d'erreur 
residuel quasi nul n'est plus vraie lorsque les caracteristiques des canaux varient dans le 
temps (canaux non stationnaires), notamment dans les reseaux sans fil et mobiles, etpour 
une complexity du code de canal realiste. En outre, l'ajout de redondance par codes 

2 5 correcteurs d'erreurs conduit a r6duire le debit utile. 

Le besoin se fait done sentir de nouvelles solutions de codage/d6codage. 
La pr^sente invention vient proposer des avanc£es en ce sens. 
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L'invention concerne un proced6 de decodage comprenant les etapes suivantes: 

a. prevoir un modele source propre a d6finir des probabilitds pour des symboles x-aires, 

associes a des transitions entre etats sources dans un diagramme d'etats source, 



un modele dit "produit", propre a definir une correspondance entre les symboles 
x-aires et de premiers elements binaires, en fonction d'6tats produits dans un diagramme 
d'etats "produit", un etat produit 6tant fonction d'un etat source 

b. recevoir des proprietes de "canal", propres a d6finir des probabilites de recevoir de 
5 seconds elements binaires, sur emission de premiers 616ments binaires, 

c. pour un flux code recu de seconds 61ements binaires, calculer des probabilites de decoder 
des symboles x-aires, connaissant des bits du flux de seconds elements binaires, 
respectivement des probabilites d' aboutir dans des etats produits connaissant desbits du flux 
de seconds elements binaires, 

10 a partir des proprietes de canal, du modele source et du modele produit, 
d. reconstimer le flux de symboles x-a^ 
a l'etape c. 

Avantageusement, le modele produit est etabli a partir du modele source et d'un modele de 
15 transfert, propre a etablir une correspondance entre des symboles x-aires et de premiers 
elements binaires, en fonction de probabilit6s associees a des symboles x-aires. Le modele 
source comprend des transitions associees a une correspondance entre des symboles m-aires 
sources et des symboles x-aires cibles, en particulier mkx et x=2. 

2 0 De preference, les modeles de produit et de transfert sont de type arithmetique et le modele 
de transfert est un modele de codage ou de d6codage. 

Plus precisement, l'etape c. comprend le calcul, par 6tapes, des probabilites de decoder un 
symbole x-aire connaissant des bits du flux de seconds dements binaires. 

Dans une premiere realisation selon l'invention, l'etape c. comprend, pour chaque etape a 
partir d'un nombre donne d'etapes et selon.un critere de selection donne, la s61ection de 
certaines probabilit6s parmi les probabilites de decoder un symbole x-aire, connaissant des 
bits du flux de seconds elements binaires, calculees lors d'une etape. 

Avantageusement, le critere de selection de l'etape c. comprend un nombre fixe de 
probabUit6s a garder parmi les plus grandes des probabilites. Alternativement ou 
additionnellement, le critere de selection de l'etape c. comprend un seuil minimal compar6 
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aux probabilites de recevoir de seconds elements binaires, sur emission de premiers 
elements binaires, de maniere a ne memoriser, & chaque etape, que les probability de 
decoder un symbole x-aire, connaissant des bits du flux de seconds elements binaires, 
calculees a partir des probabilites de recevoir de seconds elements binaires, sur emission de 
5 premiers elements binaires, plus grandes que ce seuil minimal. 

De preference, P etape c. comprend F elaboration d'un arbre constitue d'etats reliees par des 
transitions, chaque etat correspondant a une etape de calcul de probabilites et chaque 
transition, au depart d'un 6tat, correspondant a Tune de ces probabilites calculees. Ainsi, la 
1 0 selection de certaines probabilites parmi les probabilites de decoder un symbole x-aire, 
connaissant des bits du flux de seconds elements binaires, correspond a la construction 
d a 6tats correspondant & un sous-ensemble de transitions de Farbre fonction desdites certaines 
probabilites. 

i. 

15 De plus, P etape d. comprend le calcul du produit des probabilites correspondant & des 
transitions successives dans Parbre et la selection du produit le plus grand correspondant aux 
transitions successives les plus probables et au flux de symboles x-aires le plus probable; 

Selon une premifere option, F6tape c. comprend le calcul, pour un nombre donne de 
20 symboles x-aires successifs a decoder, des probabilites de decoder une sequence de 
symboles x-aire connaissant des bits du flux de seconds Elements binaires. 

Selon une deuxieme option, F etape c. comprend le calcul, pour un nombre donn£ de bits 
successifs du flux de seconds Elements binaires, des probabilites de decoder une sequence 
25 de symboles x-aire, connaissant des bits du flux de seconds elements binaires. 

Dans une deuxieme realisation selon F invention, k partir du modele produit etablissant ladite 
correspondance, Fetape c. comprend 

cl. la construction successive d'6tats relies par des transitions d'un diagramme d'6tats 
3 0 produit depuis un etat d'origine a un etat final, chaque 6tat definissant un nombre attendu 
de symboles x-aires decodes pour un nombre donne de bits du flux de seconds elements 
binaires depuis Fetat d'origine. 
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Avantageusement, la construction de l'etape c. comprend 

c2.1a suppression successive d'etats du diagrammed' 6tats produit, depuis 1' etat final a 1' etat 
d'origine, pour lesquels le nombre attendu de symboles x-aires pour un nombre donn6 de 
bit du flux de seconds elements binaires est different d'un nombre connu maximal de 
symboles x-aires. 

Selon une premiere option, l'etape c. comprend 

cl. le calcul, apres chaque 6tat construit, des premieres probabilites successives d'aboutir 
dans un etat donne connaissant des bits passes du flux de seconds el6ments binaires, 
c2. le calcul, pour chaque Stat depuis l'etat final a 1' etat d'origine, des secondes probabilites 
successives d'aboutir dans cet etat donne connaissant des prochains bits du flux de seconds 
elements binaires. 

Selon une deuxieme option, l'6tape c. comprend 

c3. le calcul, pour chaque etat du diagramme depuis l'etat d'origine a l'etat final, des 
premieres probabilit6s successives d'aboutir dans cet etat donne connaissant des bits passes 
du flux de seconds elements binaires, et , pour chaque 6tat depuis l'etat final a l'etat 
d'origine, des secondes probabilit6s successives d'aboutir dans un etat donn6 connaissant 
des prochains bits du flux de seconds elements binaires. 

De preference, l'6tape c. comprend le calcul de la probabilit6 de chaque etat en effectuant 
le produit des premieres et deuxiemes probabilites pour chaque 6tat et l'6tape d. comprend 
le calcul des produits de probabilites des differents etats successifs possibles puis la 
selection du produit maximal. 

Avantageusement, l'6tape d. comprend l'etablissement du flux de symboles x-aires 
correspondant a ces etats successifs en utilisant la correspondance du modele produit. 

L'mventionconceme6galementundecodeurc»mprenantunepremiereentr6epour recevoir 
0 un flux cod6 de seconds elements binaires. Ce decodeur comprenant une deuxieme entree 
pour recevoir des propri6tes de canal propres a definir des probabifites de recevoir de 
seconds elements binaires, sur emission de premiers 61ements binaires, un premier module 
definissant 
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- un modele source propre a definir des probabilites pour des symboles x-aires associes a des 
transitions entre etats sources dans un diagramme d'etats source, 

- un modele dit "produit", propre a definir une correspondance entre les symboles x-aires 
et de premiers elements binaires, en fonction d'etats produits dans un diagramme d'etats 

5 "produit", un etat produit etant fonction d'un 6tat source, 

- un module de traitement adapte pour calculer des probabilites de decoder des symboles 
x-aires, connaissant des bits du flux code de seconds elements binaires, respectivement des 
probability d'aboutir dans des 6tats produits connaissant des bits du flux de seconds ele- 
ments binaires, & partir des propriety du canal, du modele source et du modele produit, 

10 et adapte pour reconstituer le flux de symboles x-aires le plus probable a partir de ces 
probabilites calculees. 

D'autres caracteristiques et avantages de Invention apparaitront a Fexamen de la description 
d^taillee ci-apres, et des dessins annexes, sur lesquels : 

15 

- la figure 1 illustre un ensemble comprenant un codeur et un decodeur selon Finvention;. 

- la figure 2 est une vue gen6rale du proc6de de codage implementable dans un codeur selon 
Finvention, 

20 - la figure 3 est un diagramme formant vue generate du procede de decodage implementable 
dans un decodeur selon Finvention, 

- les figures 4A, 4B et 4C illustrent trois representations differentes d'un mSme modele dit 
"source" etablissant une correspondance entre les symboles de la source et des sequences 
de bits, 

2 5 - les figures 5A et 5B illustrent deux representations differentes d 'une meme correspondance 
entre des sequences de symboles et des sequences de bits, basees sur le module dit "source" 
des figures 4A et 4B, 

- la table 1 illustre une representation sous forme de table d'un modele dit "codeur", 

- la table 2 illustre une representation sous forme de table d'un modele dit "decodeur", 
30 - la figure 6A illustre une representation sous forme de table d'un modele dit "produit 

source-codeur"selon Finvention, 

- la figure 6B illustre une representation sous forme de treillis du modele dit "produit source- 
codeur" de la figure 6A selon Finvention, 




- les figures 7A et 7B illustrent deux representations partielles sous forme de trefflis du 
modele dit "codeur" de la table 1, 

- la figure 8 illustre une representation generalisee d'un modele dit "produit source- 

codeur"selon 1' invention, 
5 - la figure 9A illustre une representation sous forme de table d'un modele dit "produit 

source-decodeur"selon P invention, 

-lafigure9Bffiustremere P r6sentationsous forme detreillisdu modele dit "produitsource- 
d6codeur" de la figure 9A selon l'invention, 

- la figure 10 illustre une representation gen6ralisee d'un modele dit "produit source- 

10 decodeur"selon rinvention, 

-lafigurelimustreleprocedededecodagedansunepremierereaUsationselonl'mventi 

-lafigurel2illustreleprocedede decodage dans une deuxieme realisation selon 1' invention, 

- la figure 13 illustre une representation sous forme de treillis realisee a partir du procede de 

decodage de la figure 11 ou 12 , 

- la figure 14 fflustre le procede de de^^ 

- la figure 15 illustre le proc6d6 de decodage dans une quatrieme rSalisation selon 
rinvention. 
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En outre: 

- l'annexe 1 contient des expressions mathematiques utilises dans la presente description. 

Les dessins et les annexes a la description comprennent, pour l'essentiel, des 61ements de 
caractere certain. lis pourront done non seulement servir a mieux faire comprendre la 
description, mais aussi contribuer a la definition de l'invention, le cas ech6ant. 
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De maniere g6nerale, un systeme de compression de signaux multimedia (image, video, 
audio, parole) fait appel a des codes statistiques aussi appeles codes a longueurs variables. 
Ceux-ci permettent d>obtenir des d6bits approchant ce que les specialistes appellent 
"l'entropie du signal". Les codes les plus utilis6s dans les systemes existants (en particulier 
3 0 dans les standards) sont les codes de Huffman qui ont 6t6 d6crits dans l'ouvrage suivant : 
D.A. Huffman : "Amethod for the construction of minimum redundancy codes", Proc. IRE, 
40(1951),p.l098-1101. 
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Plus recemment, onavuun regain cTinteret pour les codes arithmetiques en reason de leurs 
performances accrues en terme de compression. lis permettent en effet de decoupler le 
processus d'encodage du modele suppose de la source. Ceci permet aisement dHitiliser des 
modeles statistiques d'ordre sup6rieur. Ces codes arithmetiques ont ete decrits dans des 
5 ouvrages tels que 

[1]- JJ. Rissanen "Generalized kraft inequality and arithmetic", IBM J.Res. Develop., 
20:198-203, Mai 1976 

[2]- JJ- Rissanen, "Arithmetic Coding as number representations", Acta Polytech. Scand. 
Math., 31:44-51, DScembre 1979 
1 0 ainsi que dans des brevets americains US 4 286 256, US 4 467 317, US 4 652 856, 

Jusque recemment, la conception des syst&mes de compression se faisait en supposant une 
qualite de service transport garantie. On supposait en effet que les couches inferieures du 
modele OSI incorporent des codes correcteurs d'erreurs garantissant un taux d'erreur residtiel 
15 vu de Implication quasi-nul. 

Les codes k longueurs variables pouvaient done Stre Iargement utilises malgre leur forte 
sensibility au bruit de transmission. Cette hypothese de taux d'erreur residuel quasi nul n'est 
plus vraie dans les reseaux sans fil et mobiles, dont les caracteristiques des canaux varieht 

2 0 dans le temps (canaux non stationnaires). Ce taux d r erreur r£siduel vu par le decodeur des 

signaux de source est souvent loin d'etre negligeable. Toute erreur dans le train binaire peut 
engendrer une de-synchronisation du decodeur et done une propagation des erreurs sur la 
suite des informations decodees. 

25 Pour pallier ce probleme de propagation, les standards de premieres generations (H.261, 
H.263, MPEG-1, MPEG-2) ont incorpore dans la syntaxe du train binaire transmis des 
marqueurs de synchronisation. Ce sont des mots de code longs (16 ou 22 bits constitues 
dWe suite de 15 ou 21 bits a '1' suivis dHin X) 1 ) non emulables par des erreurs se produisant 
sur les autres mots de code et qui peuvent done etre reconnus par le decodeur avec une 

3 0 probability proche de T. 

Cela conduit a structurer le train binaire en paquets dfilimites par ces marqueurs de 
synchronisation. Cela permet de confiner la propagation des erreurs au sein du paquet. 



1 er depot 
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Cependant si une erreur intervient en debut de paquet la suite du paquet peut etre perdue. En 
outre, la periodicite de ces marqueurs de synchronisation doit etre restreinte pour eviter une 
perte trop grande en efficacite de compression. 



Les nouveaux standards ( H.263+ et MPEG-4) ont alors fait appel a des codes a longueurs 
variables r6versibles (RVLQ. La particularite de ces codes est qu'ils peuvent etre decodes 
du premier vers le dernier bit d\in paquet, et, a Tinverse, du dernier vers le premier bit du 
paquet. 

Si une erreur s'est produite en milieu de paquet, cette symetrie du code permet de confiner 
la propagation des erreurs sur un segment en milieu de paquet au lieu de la propager jusqu'a 
la fin du paquet delimite par un marqueur de synchronisation, Cependant, la symetrie du 
code engendre une perte en efficacite de compression par rapport a un code de Huffman de 
Tordre de 10 %. En outre les codes k longueurs variables r6versibles n'evitent pas 
completement le probleme de propagation des erreurs : si une erreur se produit en debut et 
en fin de paquet, tout le paquet risque d'etre errone. 

La conception de proc6des de codage/decodage robustes au bruit de transmission et bases 
sur une famille de codes performants en compression (i.e., qui pennettent d'approcher 
Tentropie de la source), constitue un enjeu important, notamment pour les futurs systemes 
de communication multimedia (image, video, audio, parole) mobiles. Pour ces systemes, de 
nouveaux standards sont k 1 'etude a la fois au sein de 1TTU (International Telecommunica- 
tion Union) et de 11SO (International Standard Organization). 

Bien que les standards occupent une place prepond6rante dans le secteur des telecommuni- 
cations, une telle famille de codes peut aussi trouver des applications sur des marches de 
niche faisant appel a des solutions proprietaires. 

U invention d6crite conceme notamment un dispositif codeur/d6codeur pom une famille de 
codes largement utilises dans les systemes de compression comme des codes de type 
arithmetique comprenant les codes arithmetiques et les codes quasi-arithmetiques. 
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Les principes du codage arithmetique et quasi-arithmetique sont developpes ci-apres. lis sont 
developpes plus en detail dans Pouvrage [2] et les ouvrages suivants 
[3] - L H. Witten and R. M. Neal and J, G. Cleary . Arithmetic coding for data compression. 
Communications of the ACM, 30(6):520-540, Juinl987, 
5 [4] P. G. Howard and J. S. Vitter, Image and Text Compression, pages 85-112, Kluwer 
Academic Publisher, 1992. 

Dans le reste de la description, Futilisation du mot "sequence" designe un flux ou une suite 
d'elements. Par exemple, une sequence de bits, respectivement de symboles designera un 
1 0 flux ou une suite de bits, respectivement de symboles. 

Le principe du codage consiste de fagon generate a procfider aux Stapes suivantes : 

soit A = A a ... A u une sequence de symboles de source quantifies qui preiment leur valeur 

dans un alphabet fini A compose de M = 2 q symboles, A={a t ... a M }. Cet alphabet a une 

15 distribution stationnaire . qui indique pour chaque symbole une probability stationnake 
associee. Cette sequence de symboles M-aires (appelee egalement source A) est convertie 
en une sequence de symboles x-aires (avec x^M). Dans le cas ou cette s6quence de 
symboles M-aires est convertie en une s6quence de symboles binaires 8=8!. ..S K , (appelee 
egalement source S), K est definit par K=qxL. Les variables M et x definissent des enters 

2 0 a partir de 2. Ainsi, dans Pexemple, x = 2 et M=4. 

A titre d'exemple uniquement dans le reste de la description, et pour des raisons de 
simplification, on considerera que x=2. Cette conversion utilise un modele de source decrit 
ci-apres comprenant une correspondance entre la source M-aire A et la source binaire 5. 

2 5 Cette sequence de symboles binaires S est convertie a son tour en une sequence de premiers 

Elements binaires dHnformation JJ-V t ...U N , par lintermediaire d\m codeur qua- 
si-arithmetique de precision T. La longueur N de ce train binaire est une variable aleatoire, 
fonction de S, la sequence de symboles binaires, elle-mSme fonction de A, la sequence de 
symboles M-aires. Le train de premiers etements binaires U est transmis, a titre d'exemple 

3 0 uniquement, sur un canal sans memoire. Ce train de premiers elements binaires U est regu 

comme un train de seconds elements binaires Y en raison d'erreurs introduites lors de la 
transmission par le canal de communication. Les couches infSrieures du sysfeme de 
communication permettent en effet de determiner au moins les proprietes courantes Ic du 




canal de communication, par exemple le taux d'erreurs introduces par le canal dans les bits 
communiques, des mesures de confiance sur des bits regus. Ces dernieres mesures 
pennettent d'estimer des probabilites de recevoir un second element binaire connaissant un 
premier element binaire communiqu6 sur le canal, Ces probabilites sont estimees comme 
indique en annexe (equation 21 avec Pindice i etant un entier) a partir d'une des proprietes 
courantes Ic du canal : p, la probability d'erreur sur le canal. 

Le principe du decodage consiste k estimer la suite de symboles A, connaissant les seconds 
elements binaires appeles egalement par la suite "observations y" Dans le reste de la 
description, les lettres majuscules seront utilisees pour denoter les variables aleatoires, et les 
minuscules pour les realisations, ou observations de ces variables* Pour noter une suite de 
variables successives, la notation X u v X^...?^}. 

Le principe de codage/dfcodage peut egalement etre directement applique & une source 
correspondant non pas a des symboles M-aires (A) avec M>2 mais k des symboles binaires 
(5) avec M=2. Dans la suite de la description, un symbole x-aire, respectivement un symbole 
M-aire, comprend notamment le cas du symbole binaire (x=2, respectivement M=2) et le cas 
du symbole pour lequel x>2, respectivement M>2. 

La figure 1 illustre un exemple d'ensemble de codeur 1 et de d6codeur 2. Le codeur et le 
decodeur ont respectivement une entree Ec, Ed et une sortie Wc, Wd> Lors du codage, un 
flux de symboles x-aires k coder est regu en entree du codeur en Ec. Le codeur transmet en 
sortie Wc un flux de premiers elements binaires, appelS premier flux de bits Uc, gener£ par 
codage du flux de symboles x-aires, par exemple sous forme comprimee. Ce premier flux 
de bits Uc est envoye, par P intermediate d'un canal 11, k Pentr6e Ed du decodeur 2. Le 
canal, ay ant pu erroner le premier flux de bit (par exemple inverser certains bits du premier 
flux), transmet un flux de seconds 616ments binaires, appele second flux de bits Yd, a 
P entree du decodeur 2. Selon Pinvention, un module de traitement 6 regoit et traite ce 
second flux de bits ainsi que les proprietes Ic du canal, communiques par le canal, pour 
obtenir en sortie Wd du decodeur, un flux de symboles Sd. Un module de traitement 6 est 
propre a travailler avec le module de modeles 12. Ce module de modeles comprenant 
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-un modele de source 8 permettant de definir des probabilites pour des symboles x-aires 
associes & des transitions entre etats sources dans un diagramme d'etats source decrit ci-apres, 
- un modele dit "produif 9 7 permettant de definir une correspondance entre les symboles x- 
aires et de premiers elements binaires, en fonction d'etats dans un diagramme d'etats 
5 "produit", un etat produit etant fonction d'un etat source et d'un etat d'un modele de 
transfert. 

Le modele produit 7 est etabli a partir du modele source 8 et d'un modele de transfert 9 ou 
10. Le modele produit peut etre construit dynamiquement a partir du modele de source et des 
1 0 modeles de transfert 9 ou 10. Dans une autre realisation possible, le modele produit peut etre 
impl6mente dans le module de modeles sans etre construit dynamiquement. Le modele de 
transfert comprend un modele de codage 9 ou un modele de decodage 10. Ce modele de 
transfert permet d'etablir une correspondance entre des symboles binaires et de premiers 
elements binaires, en fonction de probabilites associees a des symboles binaires. 

15 

II est egalement possible de partir d'une sequence de symboles M-aire pouvant etre au 
pr6alable transformer en une sequence de symboles x-aires avant d'etre compress6e pn 
utilisant le modele source. 

2 0 Dans une realisation possible, ces modeles peuvent etre pre-calcul€s et mis en memoire dans 
un module externe au module de modele du decodeur. lis peuvent egalement etre utilises k 
distance, par exemple en se servant d'un reseau local de type Ethernet ou un r6seau plus 
6tendu de type internet. 

25 Le module de traitement 6 est adapte pour calculer des probabilites de symboles x-aires, 
connaissant les bits du flux de seconds elements binaires, k partir des proprietes Ic du canal, 
du modele source et du modele produit. Le module de traitement 6 est adapte pour 
reconstituer le flux de symboles x-aires leplus probable a partir de ces probabilit6s calculees 
et pour reconstituer le flux de symboles M-aires a partir du modele source si un flux de 

30 symboles M-aires a ete cod£. ... 
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Le precede decoulant de ce dispositif selon 1'invention est decrit dans le reste de la 
description et permettra, notaxmnent, de mieux comprendre le fonctionnement d'un tel 
dispositif. 

5 La figure 2 illustre un exemple de precede de codage pouvant utiliser 1'insertion de 
marqueurs de synchronisation. C^procede est implemente dans m 
une source (M-aire), par exemple une sequence de donnees. 

A I'etape 100, lorsque la sequence de donnees ne n6cessite pas de marqueurs M-aires, le 
1 o precede continue a V etape 104. Si, au contraire, des marqueurs M-aires doivent etre inser6s 
dans la sequence de donnees, la frequence, la taille et la valeur de ces marqueurs sont 
fournies par le codeur pour effectuer l'insertion de ces marqueurs a I'etape 102. 

A 1' etape 104, il est verifi6 si la sequence de donnees (comprenant ou non des marqueurs de 
15 S ynchronisation)n6cessiteuneconver S ionM-ake/binaire.Siteln'estpas 

de donnees est done binaire et est soumise au codage arithm6tique (quasi-arithm6tique dans 
rexemple de la figure) a I'etape 112. Sinon, la conversion M-aire/binaiie de la sequence 
s'effectue a I'etape 106. A I'etape 108, si des marqueurs binaires doivent etre inser6s, la 
£r6quence, la taille et la valeur de ces marqueurs sont fournies par le codeur pour effectuer 
2 0 l'insertion de ces marqueurs a l'6tape 110. 

A l'6tape 112, la sequence de donnees M-aire ou binaure est soumise au codage quasi- 
arithmetique dont la complexit6 est control par le parametre T. Si T est grand, le codage 
quasi-arithmetiquedevie*^^ 
2 5 maximale de la sequence de donn6es binaire en un train de bits. 

La figure 3 illustre un exemple de precede de decodage prenant en compte les marqueurs de 
synchronisation ins6r6s. Ce procede est implement dans un d6codeur recevant en entree le 
train binaire et des informations concernant la source M-aire ou x-aire (notamment binaire), 
30 la fr6quence, la taille et la valeur des marqueurs de synchronisation eventuels. 



A I'etape 120, differents modes de d6codage peuvent 6tre selectionnes : 
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- un decodage instantane i Fetape 124, 

- un decodage sequentiel selon F invention, qui est detaille dans le reste de la description, a 
Fetape 122, 

- un decodage optimal, qui est detaille dans le reste de la description, selon Finvention a 
5 l'Stape 126. 

A Fetape 128, il est verifi6 si la sequence de donnSes obtenue apres decodage necessite une 
conversion binaire/M-aire. 

10 Si tel est le cas, la sequence de donnees est done binaire et a Fetape 130, il est verifie si la 
sequence de donnees comprend des marqueurs binaires. Si tel est le cas, ces marqueurs sont 
supprimes k Fetape 132 en utilisant les informations sur leur frequence et leur taille. Apres 
Fetape 130 ou Fetape 132, la sequence de donnees est soumise a la conversion binaire/M- 
aire a Fetape 134. 

15 < 
Apres Fetape 128 ou 134, il est verifie si la sequence de donnees comprend des marqueurs 
M-aires a Fetape 136. Si tel n'est pas le cas, le decodage de la sequence de donnees (<>u 
source) est termine, sinon, les marqueurs sont supprimes a Fetape 138 pour obtenir la 
sequence de donnees (ou source). \ 

20 

Les figures 4 a 5 illustrent des modeles de source pour la conversion M-aire/binaire et 
binaire/M-aire d'une sequence de donn6es (ou source). 

La figure 4A illustre la representation binaire des 4 symboles d'une source 4-aire, c'est-li- 

2 5 dire une source prenant valeur dans un alphabet 4-aire. Ce modele source peut etre un arbre 

binaire de profondeur q, comme cela est montr€ sur la figure 4A pour q=2 et M=4. Dans une 
variante, cet arbre peut etre reprgsente par xme table de correspondance. Le noeud noir est 
appele racine de Farbre, les noeuds blancs sont les noeuds interm6diaires, les noeuds grises 
sont appeles feuilles de Farbre. Les probabilites P de chaque symbole M-aire correspondent 

3 0 aux probabilit6s stationnaires de la source. Les probabilites de transition (non representees) 

sur les branches de cet arbre binaire sont calculees a partir des probabilites stationnaires. 
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Ufigure4Brepresent^^ 

independants en une s6quence de symboles binaires. Ce modele de source pent egalement 
6 trere P resent6enfigure4B^ 

initial C=l et C=2 les 6tats mtermediaires. En variante, ce modele de source pent etre 
represent en figure 4C sous forme d'un treillis avecles memes etats : C=0 a gauche, 1> 6tat 
initial relie aux etats intermedials Ol et C=2 a droite, ces 6tats intermediaires, a gauche, 
etant relies a 1'etat 0=0 a droite. De plus, dans cet exemple, les probability de chaque 
symbole x-aire (ici, symbole binaire) sont indiquees sur les transitions correspondantes. 

Pour un flux de donnees comprenant plusieurs symboles M-aires non ind6pendants, un 
modele complet de la source est obtenu en connectant des modeles locaux successifs 
comme indique sur la figure 5A. Un modele local correspond a un modele de source pour 
un flux comprenant un symbole comme illustre sur la figure 4A. En reference aux figures 
5 pour un flux de deux symboles M-aires dependants, le modele complet de la source 
15 consiste a identifier les feuilles du premier arbre binaire avec la racine de 1'arbre suivant, 
racine representee par un noeud noir sur les figures 5A et 5B. La figure 5A combine amsx 
des arbres binaires entre-eux. En variante, la figure 5B combine un premier arbre binanre 
avec des automates corresponds prenant pour 6tat initial et final les feuilles de 1'arbre 
binaire. Cet automate permet de convertir une sequence de symboles M-aires dependants en 
2 0 une s6quence de symboles binaires. 

D'apres la figure 5A, les 6tats C, d6finissent les noeuds .possibles dans lesquels peut etre 
1'automate apres la production de k symboles binaires. Cet ensemble de noeud possible est 
' symbolise sur la figure par un trait mixte reliant les noeuds d'un meme etat k. La sequence 
25 Co,...,^ est une chaine de Markov etle flux de symbolesbinaires correspondant est fonction 
del transitions de cette chaine. Les transitions appelees S k puisqu'elles engendrent 
remission des symboles binaires S k , represented V ensemble des transitions possa>les entre 
l'ensemble des etats Q, a et l'ensemble des etats Q du modele source. 

30 A partir de cet automate stochastique, le flux de symboles binaires peut etre modelis6 par 
un modele de Markov cache. 
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Connaissant chaque probabilite stationnaire associee a chaque symbolebinaire, chacune de 
ces probabilites correspondant a un des segments de l'intervalleunitaire [0,1 [, le codage 
arithmetique permet, pour une sequence de symboles binaiies, de detenniner 1* intervalle de 
cette sequence. Ce dernier represente un segment de 1' intervalle unitaire et sa borne 
5 inf erieure est issue de la somme des probabilites des sequences de symboles correspondant 
aux sous-intervalles inferieurs. La largeur de cet intervalle est donnee par la probabilite de 
la sequence de symboles correspondant a cet intervalle. On appelle "intervalle courant", 
I'intervalle correspondant a une partie de la sequence et done a Tetat courant de cette 
partie de sequence, lors du calcul de V intervalle ou sous-segment associe a la sous sequence 
10 correspondant k Fetat Q. 

Des mises en oeuvre pratiques du codage arithmetique ont d'abord et6 introduites dans 
Pouvrage [1], l'ouvrage suivant : 

[5]- R. Pasco. Source coding algorithms for fast data compression. PhD thesis, Dept. of 

15 Electrical Engineering, Stanford Univ., Stanford Calif. ,1976 

et revisit6es dans [2]. Un probleme qui se pose lors de cette mise en oeuvre est la n6cessite 
de disposer d'line grande precision num6rique pour repr6senter des intervalles avec des 
nombres reels tres petits. Cette difficulte peut etre resolue en se basant sur la representation 
binaire des nombres reels dans l'intervalle [0,1 [, ce qui est decrit dans [3] .Tout nombre 

2 0 appartenant a Tintervalle [0,0.5 [ a son premier bit €gal a 0,alors que tout nombre appartenant 
a l'intervalle I0.5,l[ a son premier bit 6gal & 1. Par consequent, pendant le processus de 
codage, des que l'intervalle courant est entierement contenu dans [0,0.5[ ou [0.5,1[, le bit 
correspondant est 6mis et la taille de l'intervalle est doublee. Un traitement sp6cifique est 
necessaire pour les intervalles qui chevauchent Vi . Si l'intervalle courant chevauche Vz et est 

2 5 entierement contenu dans [0.25,0.75 [, il ne peut pas etre identifie par un bit unique. Sa taille 
est tout de mgme doublee, sans emettre de bit, et le nombre de fois (n) oii cette operation est 
effectu6e avant d'emettre un bit est m6morisee. n est le nombre de remises a P6chelle 
effectuees depuis la demiere emission d'un bit. Lorsqu'un intervalle pour lequel un bit U f = 
Uj est 6mis est atteint, alors ce bit est suivi par la sequence de n fois le bit de valeur opposee. 

30 
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^utilisation de cette technique garantit que I'intervalle courant satisfait en permanence low 
< 0.25 < 0.5 < high ou low < 0.5 < 0.75 * high, avec low et high les bornes inferieures et 
sup6rieures de I'intervalle coutant respectivement. 

5 Le principe du codage arithm6tique pent se modeliser par un automate stochastique dont les 
etats E sont definis par trois variables : low, up et n^ le nombre de remises a l'echelle. 
Linconvenient de ce codage est qu'il pent exister un tres grand nombre d'etats si la 
distribution stationnaire de la source n'est pas connue a priori ou si le flux comprend un 
nombre de symboles x-aires important. Le codage quasi-arithmetique, egalement appel6 
10 codage arithmStique a pr6cision reduite, permet de r6duire le nombre d'6tats possibles sans 
perte importante en efficacite de compression, comme d6crit dans Pouvrage [4]. Pour 
effectuer le codage qnasi-arithmetique, I'intervalle r6el [0,1[ est remplace par I'intervalle 
entier [0,T[. L'entier T permet de controler le compromis entre la complexity et l'efficacitS 
de compression : si T est grand, alors les subdivisions de I'intervalle repr6senteront 
15 P r6cis6ment la distribution de la source. Par centre, si t est petit (proche de 1), toutes les 
subdivisions possibles de I'intervalle pourront 6tre pre-calculees, ces subdivisions 
represented seulement une approximation de la distribution de la source. 

Comme pr6sente dans la table 1, toutes les transitions entre 6tats et les emissions de bits 
20 correspondantes peuvent etre pre-calculees. Les operations arithm6tiques sont alors 
remplac6es par des lectures de table. 

Une sequence S* de symboles binaires est convertie un flux de bits U, N par un arbre binaire. 
Cet arbre pent etre vu comme un automate ou un modele codeur qui modelise la distribution 
25 du train binaire. Le codage dNm symbole determine le choix d'un arc de l'arbre, auquel pent 
etre associee remission de bits. Chaque noeud correspond a un etat E du codeur arithme- 
tique. Les transitions successives entre ces 6tats suivent la distribution de la source P(S k | S k _ 
,) pour une source de Markov d'ordre 1. Soit E, V6tat de l'automate a chaque "instant 
symbole" k. Comme pour le codage arithmetique, l'etat E, du codeur quasi-arithm6tique est 
30 defini par trois variables: lowS k , upS k et nscl k . Les termes lowS k et upS k d6signent les 
bornes de I'intervalle resultant des subdivisions successives de [0,T[ correspondant au 
codage de la s6quence S a K . La quantit6 nscl k est remise a zero chaque fois quHm bit est emis 
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et incrementee chaque fois qu'une remise a Techelle sans emission de bit est effectuee. Elle 
repnSsente done le nombre de remises & Techelle effectuees depuis la demiere emission de 
bit. Quand un bit est emis, il est suivi par nscl k bits de valetir opposee. 

5 Dans la mesure ou il y a un nombre fini de subdivisions possibles de l'intervalle [0,T[, il doit 
etre possible de pre-calculer tous les etats du codeur quasi-arithmetique sans connaitre la 
source a priori. 

A titre d'exemple, la table 1 prSsente les 6tats, les sorties et toutes les transitions possibles 
1 0 dHin codeur quasi-arithmetique de precision T=4 pour une source binaire. On s'interessera 
dans un premier temps aux ensembles CI et C2 de colonnes. La valeur de la variable nscl 
k n'est pas inclue dans le modele d'etat Seuls les increments de cette variable sont signales 
par la lettre f dans le tableau comme dans Fouvrage [4]„ 

15 Le codeur a trois etats possibles, correspondant a des subdivisions entires de rintervalld. 
[0,4[ prises par l'iiitervalle courant [lowS k , upS k [ . On se reportera a I'ouvrage \ 

[6] P. G. Howard and J. S. Vitter, Design and Analysis of fast text compression based on 
quasi-arithmetic coding. In Data Compression Conference, pages 98-107, Snowbirdj 
20 Utah,Mars-Avril 1993. 

pour le choix de ces subdivisions. En fonction de la probabilite du symbole binaire d'entree 
0, l'intervalle courant [lowS k , upS k [ sera subdivis6 en intervalles entiers. Suivant que, dans 
la colonne C2, le symbole binaire d?entree est 0 ou 1, les subdivisions en intervalles entiers 
25 de l'intervalle courant PowS k , upS k [ possibles conduisent a remission de bits indiquees ou 
a la remise a l'echelle appropriee avant de determiner le prochain etat de 1'automate. 

Dans le cas ou les symboles en entree sont des symboles binaires, le nombre d'etats peut etre 
encore reduit en identifiant les symboles binaires comme " moins probable" (LPS) et "plus 
3 0 probable" (MPS) dans la colonne C3, plutot que 0 et 1 dans la colonne C2. Cela revient a 
considerer la probabilite de chaque symbole binaire plutot que la valeur du symbole binaire, . 
les MPS et LPS pouvant corresponds soit a un bit 0 soit a un bit 1. Le nombre de 
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combinaisons possibles des probabilites de la source binaire est reduit Cela pennet de 
combiner certaines transitions et d'eliminer des etats de P automate. 

La sequence E^..!^ d'6tats successifs du codeur quasi-arithmetique forme une chaine de 
5. Markov et les sorties du codeur sont fonctions des transitions sur cette chaine. 

De la meme maniere que le codage quasi-arithm€tique, le decodage peut se modeliser par 
un automate stochastique & nombre d'etats E n fini comme presente sur la table 2. 

10 L'etat d'un decodeur arithm6tique est d£5ni par deux intervalles : - [lowU D ,upU n [ et 
powSg^upSj^ [. L*intervalle powU n ,upU n [ correspond au segment de Tintervalle [0,1[ defini 
par la reception de la sequence de bits U t n . L'intervalle PowS^upSi^f est le segment de 
l'intervalle [0,1[ obtenu quand la sequence de symboles a ete decodee. Cet intervalle est 
mis a jour par cumul des probabilites de chaque symbole binaire decodfi, de la meme 

1 5 maniere qu'au codage pour chaque symbole code. II reproduit le f onctionnement du codeur. 

Un codeur et un decodeur arithmetique ne sont pas synchronis6s. En effet, au codage, la 
reception d'un symbole n'entraine pas n£cessairement remission d'un bit De meme, au 
d6codage, la lecture d'un bit n'entraine pas necessairement la production d'un symbole. 
2 0 Ainsi, si la lecture du premier symbole de la sequence a coder n'entraine aucune emission 
de bit, il se peut qu'au decodage, la lecture du premier bit suffise a decoder le premier 
symbole. Les deux intervalles sont mis & l'echelle comme au codage afin d'eviter les 
problemes de precision numerique. 

2 5 Dans Pexemple de la table 2, un automate stochastique est defini pour une valeur de T = 4* 

La table pr6sente les 6tats, transitions et sorties d'un decodeur quasi-arithmetique pour une 
source binaire, avec T=4 et la simplification MPS/LPS. 

Soit E n l'etat de cet automate a "1 'instant bit" n, c'est-a-dire a P instant de traitement du bit 

3 0 regu U n . Cet etat est defini par les valeurs des intervalles powU D ,upU n [ et PowS^upSj^ [. 

En fonction de la distribution stationnaire de la source (representee, colonne Kl, par la 
probability du symbole MPS le plus probable), l'intervalle courant PowS^upS^ [ sera 
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subdivise en intervalles entiers. Suivant que le bit d'entr6e est U n = 0 dans la colonne K2 ou 
U n =1 dans la colonne K3, les subdivisions en intervalles entiers de Pintervalle courant 
powS^upS^ [ possibles conduisent a remission de symboles MPS ou LPS indiqu6es ou 
a la remise a 1'echelle appropriee avant de determiner le prochain etat de P automate. 

5 

Un codeur, respectivement decodeur utilisant P automate stochastique de la table 1, 
respectivement table 2, d&narre a T etat ^ =1, respectivement E n =1, 

Differents modeles sont envisages selon P invention pour etre utilis6s dans des procedes de 
10 codage/d6codage et pour etre implementes dans les codeurs/decodeurs correspondants. 
Ainsi, il est envisage un modele dit "produif ' pouvant etre un modele resultant du produit 
des modeles source et codeur ou un modele resultant du produit des modules source et 
decodeur, 

15 De maniere g6nerale, les etats du modele produit doivent rassembler les informations 
contenues dans les etats de chacun des deux modeles (source et codeur ou source et 
decodeur). Dans le cas general, la taille de Tespace d'etats du "modele produit" ne peut p]as 
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etre connue a priori. Elle peut etre calculee a partir du modele de source binaire et 
parametre T. 
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Pour un modele produit source/codeur, comme le nombre de bits, dans chaque element 
binaire, produit par chaque transition du modele produit est aleatoire, la structure des 
dependances entre la sequence de seconds elements binaires et les 6tats du modele produit 
est aleatoire. Afin de "capturer"' cette structure de dependances, le processus de Markov, 

2 5 represents par les etats X* du codeur selon le modele produit, est augmente d'un compteur 

N k pour obtenir la structure de dgpendance suivante (X,N)=(X 1 ,N 1 )... (X^N*). Chaque etat 
X k du modele produit est defmi par Petat du modele codeur, P etat du modele source et le 
nombre de bits emis quand le k-i&me symbole a 6t6 code : (X^ISQ = (lowS^upS^Ck^STk). 
Cela conduit & la structure de dependances decrite graphiquement par la figure 8. Suivant 

3 0 cette structure, une sequence de symboles binaires S X K est convertie en une sequence de bits 

U^avec N K le nombre de bits 6mis.quand K symboles ont 6te cod£s. Etant donn£s un etat 
X M et un symbole d'entree S^, le modele produit specifie les bits k U m qui doivent 
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etre emis et l'etat suivant X k . Certaines transitions peuvent n'emettre aucun bit. Les 
probabilites des transitions entre (X^NJ = OowS^upS^CW et l'etat suivant dans le 
treillis (X k+1 ,N k+1 ) = CowS k+1 ,upS k+1 A +1 ,N k+1 ) sont donnees par le modele de source binaire, 
PCQjC^ PCS^IQO. Les seconds elements binaires a Y w sont obtenues a la sortie 

5 du canal de transmission lorsque les premiers elements binaires V m . M a U m ont ete emis. 
Ces seconds elements binaires sont des observations sur les premiers el6ments binaires. 

Un exemple de modele produit issu des modeles source et codeur est presente d'apres les 
figures 6 et 7. L'exemple de modele produit definit dans la table de la figure 6A est issu du 
10 produit du modele source de la figure 4C et du modele codeur de la table 1. En theorie, la 
taille de l'espace d'&ats r6sultant de ce produit est donnee par le produit des tailles des 
espaces d'Stats des deux modeles (source et codeur). Toutefois, des simplifications ont lieu, 
ce qui permet d'aboutir a un nombre d'etats moins important. 

15 Le modele codeur aritbm6tique de la table 1 <T=4) utilisant le modele d'etat simplifie - de la 
colonne C3 a ete represent6 sous forme de treillis aux figures 7A et 7B, la figure 7A 
correspondant a une probabilite 0.63 £P(MPS), la figure 7B correspondant a une probability 
0.37 £P(MPS)<c 0.63. L'une oul'autre de ces figures peut 6tre utilisee suivant les probability 
de la source binaire. Les probabilites de la source binaire, r6sultant d'une conversion a partir 
2 0 d'une source M-aire, dependent de l'etat Q du modele de cette source binaire. Les transitions 
du modele produit source/codeur sont ainsi f onction de la distribution de la source. La table 
de la figure 6A est construite en fonction des Stats Q du modele de source de la figure 4C, 
de la probability des transitions S,,, les transitions S k =0 et S k =l de la figure 4C dtant 
modifi6es en S k =LPS et S k =MPS, et des transitions du modele de codeur repr6sente aux 
2 5 figures 7A et 7B en fonction de la probability de la source. La colonne Ml de la figure 6A 
ne contient que 4 6tats du modele produit definis dans la colonne M2 par les variables 
d'6tats du modele de source et du codeur correspondant. Le modele de la source ay ant 3 6tats 
et celui du codeur ayant 2 etats, le modele produit aura au maximum 6 etats. En pratique, 
des simplifications ont lieu comme expliquS ci-apres a l'aide des figures 7A et 7B. La figure 
30 7A montre que pour un etat E,=0 du codeur, il est possible de continuer dans un Stat E k =l 
du codeur pour une probabilite du modele de source 0,63<P(MPS). Par contre, la figure 7B 
indique qu'il n'est pas possible d'atteindre l*6tat E^l du codeur a partir de l'6tat ^=0. 



21 

Ainsi, dans le modele produit de la figure 6A, aucun etat du modele produit X k =4 et X k =5 
correspondant a des paires d'etats du modele source et codeur (E^l, Q=l) et (Ej==l, Q=2) 
n'existe pour la loi de probability telle 0,37<P(MPS) <0,63 du modele de source. Les 
colonnes M3 et M4 indiquent, pour le symbole en entr6e (MPS ou LPS), les bits emis et 
5 Petat suivant du modele produit. 

Pour un modele source/decodeur selon la figure 9A ou 9B, les etats de ce modele produit 
doivent rassembler les informations des etats des modeles de source et de dScodeur. Un etat 
^ de ce modele produit est done defini par X^OowU^upU^lowS^upS^CKj, etant le 

10 nombre de symboles decodes quand n bits ont ete regus. Le nombre de symboles produits 
par chaque transition du modele produit est al6atoire. Par cons6quent, la structure des 
dependances entre la sequence de second elements binaires et les symboles decodes est 
aleatoire. La chaine de Markov X 2 ...X N du modele produit source/decodeur est augmentee 
par une variable compteur pour donner la chaine de Markov (X 5 K)=pC l5 K 1 ) ... (X N ,K N ) 

15 ce qui aboutit a la structure de dependances illustree k la figure 10. Suivant cette structure, 
une sequence de bits U t N est convertie en une sequence de symboles S^, avec le nombre 
de symboles decode quand N bits ont ete regus. Etant donnes un etat X n et un bit d'entrle 
U n+1 , 1 'automate produit la sequence de symboles binaires S^ 1 ^ 1 , et specific Tetat suivaiit 
X n+1 . Les probability des transitions entre (X^KJ et (X^K^) dans le treillis dependent 

20 du modele de source binaire. Elles sont donnees par la formule en annexe (17). Comme pour 
le modele source/codeur, des simplifications d'Stats de ce modele sont possibles en fonction 
du modele du d£codeur et du modele de source utilise. Dans Pexemple de la figure 9 A ou 
9B, le mod&le source utilise est celui de la figure 4C et le modele decodeur, celui de la table 
2. 

25 

Les proc6d6s de d6codage robuste pour codes arithmetiques sont bases sur les mod&les de 
dependances decrits precedemment. 

Une premiere realisation du procede de d6codage robuste est illustr£e par la figure 11. 

30 

II est fait reference a l'annexe 1 pour ce qui conceme le developpement matMmatique 
associ6. 
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Le decodeur recoit en entree le train de bits a decoder. Le modele de source, la variable N 
represents le nombre de bits recus a decoder et la probabilit6 P(U n fYJ d'avoir 6mis un 
premier bit (ou premier element binaire) U n sachant le second bit (ou second element 
binaire) Y D observd sont des donn6es que le decodeur a pour entr6e, c'est a dire que, a titre 
5 d'exemple" seulement, ces donnees peuvent etre recues en entree avec le train de bits a 

decoder ou peuvent etre lues dans une memoire accessible par le decodeur. La probabilite 

P(U 0 /Y„) peut etre obtenue a partir des proprietes du canal. 

Le critere du MAP (Maximum A Posteriori) correspond a 1' estimation bayesienne optimale 
1 0 dun processus X a partir des observations disponibles Y comme indique par 1'equation (3) 
annexe 1. L'optimisation est effectuee pour toutes les "sequences" possibles. Dans cette 
premiere realisation, ce critere est applique a 1'estimation des etats caches des processus 
(X,N), qui peuvent etre considers comme des processus a horloge symbole, et des processus 
(X,K), qui peuvent etre considers comme des processus a horloge bit, etant donnees les 
15 sequences de bits observ6es. L'estimation peut etre faite soit a horloge bit, soit a horloge 
symbole. Pour des raisons de simplicite de mise en oeuvre, 1'estimation a horloge bit est 
privilegiee. 

Dans ce cas, 1'estimation de l'ensemble des 6tats caches (X 3 K><X l ,K 1 ) ... <X N ,K N ) est 
20 equivalente al'estimation dela s6quence de symboles decod6s correspondante S = S x ... S*. 
... Skn , sachant les observations Y, N a la sortie du canal. La meilleure sequence (X,K) peut 
. gtre obtenue a partir des probabilites locales sur les paires (X..KJ grace a l'equation (4). 

Le calcul de P(X n ,Kj Y) peut etre effectue a partir de l'equation de d6composition (5) ou a 
2 5 designe un facteur de normalisation. Les dependances de la chaine de Markov permettent 
un calcul recursif des deux termes de la partie droite de cette equation, base sur 1'algorithme 
BCJR decrit dans I'ouvrage 

[7] L.R. Bahl, J. Cocke, F. Jelinek, and, J.Raviv. Optimal decoding of linear codes for 
minimizing symbol error rate. IEEE trans, on transformation theory, 02:284-287, March 
30 ' 1974. 
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Ainsi, a l'etape 300, le treillis est tout d'abord construit en utilisant le modele source/deco- 
deur de la figure 9A par exemple et les probabilites correspondantes sont calculees pour 
chaque transition du treillis. Ainsi, a partir du modele produit, une passe avant permet la 
construction successive des etats relies par des transitions du diagramme d'etats produit 
5 depuis un 6tat d'origine a un 6tat final. Chaque etat definit un nombre attendu de symboles 
decodes pour un nombre donne" de bits pass6s du flux de seconds elements binaires depuis 
l'6tat d'origine. Durant cette meme etape, la passe avant de l'algorithme BCJR permet de 
calculer le terme P(X a =x n ,K n =k n |Y 1 n ) de la decomposition comme la somme sur l'ensemble 
des valeurs possibles de x a . l3 kn-i du produit de trois probabilites conditionnelles (equation 

1 0 (6): la probabilite recursive, la probability de transition entre les etats (x^, k,,.,) et (x„, kj 
domtee par le modele produit (comme indique en equation 17 et representant le produit des 
probabilites de transition P(S k |C k . 1 ) precedemment dans le modele source) et la probabilite 
d'avoir emis le bit U n (pouvant etre 0 ou 1) d6clenchant la transition entre les etats (x^, k^) 
et (x n , kj sachant le second element binaire Y n . Cette derniere probability est obtenue a 

1 5 partir des proprietes du canal. Le calcul recursif est initialise a l'6tat de depart (0,0) et permet 
de calculer Pp^^lYi 0 ) pour tout etat possible (x^kj avec n=l...N de lliorloge bit. Ainsi, 
cette passe avant permet la construction du treillis du modele produit et le calcul des 
probabilites Ppc^KjY,"), successives etpour chaque instant n, d'aboutir dans un 6tat dorute 
(X^K,,) connaissant des bits du flux de seconds etements binaires Yf. 

20 

A l'Stape 302, la passe arriere permet d'obtenir le terme P^jX^xJK^kJ de la 
decomposition pour tout etat possible (xJO a chaque instant bit n consecutivement de N 
a 1. Ce terme est calcule comme la somme du produit de trois probabilites conditionnelles 
(equation (7) ): la probabilite recursive initialisee pour tous les ^derniers 6tats" possibles 

25 (X^Kn), la probabilite de transition entre les etats (X„ + i, K„ +1 ) et (X,,, KJ donnee par le 
modele produit et la probabilite d'avoir emis le bit U n+1 declenchant la transition entre les 
etats (X^i, K„ +1 ) et (X„, sachant le second element binaire Y n+1 . Comme ptecedemment, 
cette derniere probabilite est obtenue a partir du modele de canal. Ainsi, cette passe arriere 
permet de calculer les probabilites P^" JX^K^k,,), successives et pour chaque instant 

3 0 n de N a 1, d'observer des bits du flux de seconds elements binaires connaissant un etat 
donne (X^KJ. 



1er depot 
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Lorsque le nombre de symboles K est connu, une contrainte de terminaison peut etre 
rajoutee lors de la passe arriere a Petape 302* Ainsi, tous les chemins dans le treillis qui tie 
conduisent pas au bon nombre de symboles Kj^K sont supprimes. Ainsi, la passe arriere 
comprend la suppression successive des etats du diagrannne d'6tats produit, depuis Petat 
5 final a Petat d'origine, pour lesquels le nombre attendu de symboles pour un nombre donne 
de bits du flux de seconds elements binaires est different d'un nombre connu maximal de 
symboles. Ainsi, la premiere passe de Petape 300 demarre a Tetat initial du treillis et 
construit successivement les differentes sections du treillis. Au cours de cette passe avant, 
il n f est pas possible de savoir quels sont les chemins qui ne respectent pas la contrainte de 
10 terminaison. Une fois la passe avant terminee, les etats de la derniere section (n=N) du 
treillis qui ne respectent pas la contrainte sont supprimes. La passe arrifere permet alors de 
supprimer toutes les branches qui n'aboutissent pas. 

A Petape 304, les probabilites obtenues des passes avant et arriere sont combinees sous la 
15 fonne d'un produit pour chaque valeur de n variant de 1 & N pour obtenir les probability 
P(X n ,KjY) pour chaque etat n. De cette maniere sont estimSes les probabilites liees aux 
transitions du treillis construit. 

A Petape 306, le produit de Pequation (4) est effectue sur les probabilites P^jKjYj 11 ) avec 
20 n variant de 1 a N. Ce produit est effectu6 pour les differentes combinaisons possibles de ces 
probabilites. Ainsi, la probabilite P(X,K|Y) la plus grande est retenue et les probabilites 
PCX^KjY) constituant ce produit permettent de determiner le chemin le plus vraisemblable 
dans le treillis. A partir du treillis et de la table correspondante du module produit 
source/decodeur, il est possible de determiner les symboles correspondant a ce chemin et de 
25 reconstituer une sequence de symboles S^. 

Dans une autre realisation possible, le treillis sur lequel Testimation est effectu6e peut etre 
pr€-calcule et memorise k partir du modele produit source/decodeur et de la longueur N de 
la sequence de bits observ6e. Ainsi, en reference & la figure 12, k partir du modele 
3 0 source/decodeur et des valeurs de K et N, le treillis est construit a Petape 400 en une passe 
avant et simplify par une passe arriere a Petape 402. Ce treillis est stocke en memoire pour 
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etre utilise a chaque arrivee de train de bits pour lequel le modele source/decodeur et les 
valeurs K et N sont les memes. 

Pour un train de bits regus comprenant N=6 bits, pour un nombre connu de symboles a 
trouver K=7 et pour un modele de source/decodeur illustre en figure 9A, la construction du 
treillis est illustre en figure 13. Ce treillis comporte, a chaque etat (X,,KJ, plusieurs valeurs 
couples (X^K,,). Par exemple, le couple (X^) peut prendre les valeurs (1, 1) ou (2, 1) 
selon que le bit emis U x a pour valeur 0 ou 1. 



10 A chaque arrivee de train de bits (ayant le meme modele, les memes valeurs de K et N que 
. le treillis), une passe avant et une passe arriere sont effectuees sur le train de bits aux etapes 
404 et 406, par exemple de facon parallels Les probabilitSs PCX^K^lYj 0 ) et P(Y n+1 N | X^KJ 
sont calculees comme indiqu6 aux equations (6) et (7). Le produit de ces probability est 
effectue a l'etape 408 pour chaque valeur de n et correspond a chaque transition du treillis. 

15 II est ensuite procede au calcul du produit des probability pour n variant de 1 a N a l'6tap!e 
410 comme indiquS en (4) pour obtenir des probability P(X,K|Y). fi est garde le maximum 
parmi ces probability P(X,K|Y) pour determiner le chemin le plus vraisemblable sur le 
treillis et la sequence de symboles correspondante. 

2 0 Selon l'exemple de la figure 13, le chemin le plus vraisemblable (obtenant la plus grande 
probabilite P(X,K]Y) ) est marque en gras. La suite des etats (Xo,Ko), ....(X^K,,) permet, a 
la lecture de la table 8A de determiner la sequence de symboles correspondant a ce chemin. 

D a ete present un procede de decodage optimal de codes quasi-arithmetiques. Ce procfide 
2 5 n6cessite la creation et la manipulation d'un treillis d'estimation, comme par exemple celui 
de la figure 13. La taille de ce treillis depend du nombre d'etats du modele produit 
source/decodeur (4 6tats sont possibles pour le modele produit source/decodeur de la figure 
9A par exemple) et de l'amplitude des valeurs que peut prendre la variable K„ a chaque 
instant bit n. Le nombre d'etats du modele produit source/decodeur est lui-m6me dependant 
3.0 du modele de source binaire et du parametre de precision du codage quasi-arithm6tique T. 
Pour regler le compromis entre l'efficacite de compression du codeur quasi-arithmetique et 
la complexitg du procede d'estimation, il suffit done de regler la valeur du parametre T. Une 
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valeur de T aussi petite que 4, par exemple, entraine un debit excedent raisonnable. Si 
cependant une compression maximale est requise, T doit avoir une valeur beaucoup plus 
grande et la valeur minimale de T requise depend de la distribution de la source. Dans la 
majeure partie de ces cas, la valeur de T est trop importante pour pouvoir utiliser le precede 
5 d'estimation optimal decrit ci-dessus. Lenombre d'6tats dumodele produit source/decodeur 
peut 6galeroent etre trop important pour pouvoir pre-calculer ces 6tats. 

Un procede d'estimation applicable au codage arithmetique et au codage quasi-arithmdtique 
est maintenant decrit quelle que soit la precision T choisie, e'est-a-dire quelque soit 
10 l'efficacite en compression choisie. 

Ce proced6 s'appuie sur les meme modeles de source, de codeur et de decodeur que ceux vus 
prec6denmient.I^nombred'etatsdesmodelesdecodemetdedecodeuretanttropimpo 

pourqueceux-cipuissentetreetreprS-^^ 

15 pas non plus etre construit. Le procede d'estimation propose ci-apres est done bas6 sur 
l'exploration dhin sous-ensemble du treillis d'estimation, qui est consider* comme un arbre. 
Ce procede est appel6 procedS d'estimation sequentiel. II peut s'appliquer a partir des 
modeles source/codeur et effectuer les calculs s6quentiels en utUisant "l'horloge symbole", 
e'est-a-dire un compteur incremente selon le nombre de symboles a trouver. II peut 

20 6galement s'appliquer a partir des modeles source/decodeur et effectuer les calculs 
s6quentiels en utilisant "V horloge bit", e'est-a-dire un compteur incr6mente selon le nombre 
de bits a decoder. Ce deuxieme proced6 d'estimation est pr6sente pour une source binaire a 
titre d' exemple uniquement. II peut 6galement s'appliquer a une source M-aire (M>2). 

25 La figure 14 illustre le proced6 d'estimation sequentiel dit a horloge symbole, e'est a dire 
base sur le modele source/codeur. Ainsi, d'apres la figure 8 repr6sentant la structure de 
dependance du modele source/codeur, 1' estimation de la s6quence des 6tats (X, N) pour un 
train de bits est equivalent a l'estimation de la sequence de transitions entre ces 6tats, e'est 
a dire a l'estimation de la sequence S sachant les observations Y a la sortie du canal. Ceci 

30 se traduit par une formulation mathematique present6e en (8) ou a designe un facteur de 
normalisation. Dans le cas du codage aritbm6tique optimal (pour une valeur de T elevee), 
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comme les bits emis sont distribu6s uniformement et independamment, ce facteur de 
normalisation est VQJ^/P^^) = 1. 



A l'etape 500, la probabilite de d6coder une sequence de symboles (de 1 a k) connaissant 
5 des bits du flux de seconds elements binaiies P(S f/Yj**) pour chaque instant symbole k peut 
etre calculee en une passe avant unique grace a la d6composition de liquation (9 ) ou N k est 
le nombre total de bits emis lorsqu'on est dans l'etat Xfc. En tenant compte du modele de 
source binaire defini a 1' equation(l) et en considerant un canal sans m6moire, la probabilite 
d'observer les seconds elements binaires correspondant a un etat du modele ne depend que 
10 des premiers elements binaires pour cet etat, cette probabilite peut etre reecrite comme 
indique a l'equation (10). 

Le calcul est initialise comme indique en (11) avec Sj = s x pouvant prendre deux valeurs 
possibles 0 ou 1 dans le cas d'une source binaire, et M valeurs possibles dans le cas d'une 
15 source M-aire. C x est l'etat de depart du modele de source. Pour un canal consid6r6 sans 
memoire, le second terme de cette equation modelisant le bruit introduit par ce canal 6st 
donne a l'6quation (12) comme le produit, pour un 6tat X k donnd du modele source/codeur, 
des probability d'observer un second element binaire Y n pour un premier element binaire 
transmis U n . Partant de cette initialisation, l'equation (10) est utilisee pour calcujer 

2 0 P(Si k /Y, Nk ) pour chaque instant symbole k consecutif . 

Toutes les sequences de symboles possibles 6tant considerees, la taille du treillis 
d'estimation croit exponentiellement avec le nombre de symboles K. Ainsi, l'6tape 502 
autorise le calcul des probabilit6s P(S 1 k /Y 1 Nk ) pour un nombre de symboles k variant de 1 a 
25 xl. Pour les symboles a partir de cette valeur xl, un elagage intervient a l'6tape 504 pour 
maintenir une complexity du treillis d'estimation aun niveau desire, xl estun entier inferieur 
ou egal au nombre de symboles K. 

L'elagage de l'etape 504 peut comporter differentes realisations. A titre d'exemple, deux 

3 0 r6aUsations sont presentees ci-apres. Les realisations peuvent etre combinees ou employees 

separement. 



28 



Ainsi, dans une premiere realisation de l'elagage de 1' etape 504, seuls les W noeuds les plus 
vraisemblables d'un etat (X k , N*) sont s61ectionnes et ceci a partir d'un nombre k = xl pour 
garantir la significativite des probabflit6s calculees (appel6es vraisemblances). La valeur de 
W depend des ressources de memoire disponibles. Ainsi, pour chaque etat (X k , NJ, seules 
5 lesWprobabiliteslesplusgrandespar^^^ 

de symboles (de 1 a k) connaissant des bits du flux de seconds elements binaires PCSjVY, 1 *) 
sont s61ectionnees et m6morisees. Dans une variante, il est possible de garantir une 
terminaison correcte de l'estimation en raffinant cette contrainte. Ainsi, seuls les w noeuds 
les plus vraisemblables seront selectionnes et conserves pour chaque valeur possible de N k 
1 0 (correspondant a l'horloge symbole). N k est le nombre de bits emis pour k symboles codes 
lors du codage. La relation entre W et w est donnee en (19) avec et N^ dependant 
respectivement de la probability de la s6quence S, k la plus et la moins vraisemblable. 

Dans une deuxieme realisation de l'61agage de 1' etape 504, le critere d'elagage est bas6 sur 
15 un seuil applique sur la vraisemblance des chemins. La vraisemblance d'un chemin est 
donn6e par l'equation (8) (horloge symbole). Dans cette equation, les termes P(Y|U) sont 
lies au modele de canal comme indique dans l'equation (12). P^yjU^uJ represente la 
probability qu-un bit y n soit observe sachant quW bit a 6t6 transmis ou la probability de 
recevoir un bit y n sachant q'un bit u, a 6X6 6mis, ce qui equivaut a la probability que le bit. 
2 0 observ6 soit errony. Soit p la probability d'erreur binaire moyenne qui est fixee arbitrate, 
ment comme toiyrable sur le canal durant l'estimation. L'elagage consiste a supprimer tons 
Jes chemins pour lesquels la probability d'erreur binaire moyenne a la sortie du canal est 
superieure a p. Ainsi, a chaque instant symbole N k , le critere d'yiagage comprend un seuil 
p compary aux probabilites PCY^U^ de manure a ce que tons les chemins du treillis 
25 d'estimation pour lesquels est vyrifie P^U^ < (l-p) 1 *, soient supprimys. Seules sont 
mymoriseeslesprobabmtysdedecoderunesyquencedesymboles(delak)connaissantdes 

bits du flux de seconds yHments binaires P(S 1 k /Y J Nk ), calculees a partir des probabilites 
. POT^Um,) plus grandes que le seuil minimal p. 

30 Aprcs incrementation de la valeur k d'une unite, la probability. P^/Yj 1 *) pour le nombre 
de symboles Nk (avec k>x) est calcu^e a l'elape 506. L'elagage realisy a l'etape 504 est 
rytery pour ce nouvel instant symbole Nk. Les etapes 504 et 506 sont effectuees pour une 
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valeur k incr&nentee jusqu'a la valeur K. A Petape 510, une fois la valeur K atteinte, les 
probabilites P(S 1 K /Y 1 NK ) marginales a posteriori sont connues pour toutes les s6quences s a K 
possibles. La sequence s x K la plus vraisemblable est choisie, c'est-a-dire la sequence pour 
laquelle la probabilite P(S 1 K /Y 1 NK ) est la plus grande parmi les sequences possibles du treillis 
5 d' estimation obtenu. 

Lorsque le nombre de bits N est connu, il peut constituer une contrainte sur la terminaison 
du treillis d'estimation et les chemins qui ne respectent pas la contrainte N ne sont pas 
consideres. 

10 

La figure 15 illustre le proced€ d'estimation sequentiel dit a horloge bit, c'est a dire base sur 
le modele source/d6codeur. Ainsi, d'apres la figure 10 repr6sentant la structure de 
dependance du modele source/decodeur, P estimation de la sequence des etats (X, K) pour 
obtenir un train de symboles est equivalent a P estimation de la sequence de transitions entre 
1 5 ces 6tats, c'est a dire a Testimation de la s6quence S sachant les observations Y a la sortie du 
canal. Ceci se traduit par une formulation mathematique present6e en (13) ou a designe un 
facteur de normalisation. Dans le cas du codage arithm6tique optimal (pour une valeur de 
T Slevee), les bits emis etant distribu6s uniformement et independamment, ce facteur de 
normalisation est P(U 1 °)/P(Y 1 n )=l. ; 

20 

A P6tape 600, la probabilite de decoder une s6quence de symboles connaissant les bits 
du flux de seconds elements binaires P(S 1 Kn /Y 1 n ) pour chaque bit n (defini ulterieurement) 
peut etre calculee en une passe avant unique grace a la decomposition de liquation (13 ) ou 
est le nombre total de symboles decodes a la reception du bit n dans T6tat X n . En tenant 

25 compte du modele de source binaire qui definit en (2) les probabilit6s de transition entre 
etats (X^KJ et (X^K^, la probability P(S 1 Kn ) d'obtenir une sequence de symbole 
s'exprime de maniere recurrente selon liquation (14) en utilisant Pequation (15). En 
consid£rant un canal sans memoire, c'est-a-dire que la probabilite d'observer les seconds 
61€ments binaires correspondant a un etat du modele ne depend que des premiers 61ements 

3 0 binaires transmis pour cet etat, la probability ^(Yf/Uf) peut Stre re6crite et permet d'obtenir 
la probability P(S 1 Kn /Y 1 n ) comme indiqu6 en (16). 
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Comme dans le cas de la figure 14, le calcul de la probabilite V(S«*/Y X °) par l'equation (16) 
est initialis6e pour n=l, S, = s, pouvant prendre deux valeurs possibles 0 ou 1 dans le cas 
d'une source binaire, et M valeurs possibles dans le cas d'une source M-aire. Q est l'etat de 
depart du modele de source. Partant de cette initialisation, liquation (16) est utilisee pour 
5 calculer P(S 1 Kn /Y 1 n ) pour chaque instant bit n consecutif . 

Toutes les sequences de symboles possibles etant considerees, la taille du treillis 
d'estimation croit exponentiellement avec le nombre de symboles K pouvant etre estime. 
Ainsi, l'etape 602 autorise le calcul des probabilites P(S 1 Kn /Y 1 n ) a Petape 600 pour un 
1 0 nombre de bits n variant de 1 a x2, x2 etant un entier designant un nombre seuil de bits. Pour 
les bits a partir de cette valeur x2, un elagage intervient a l'etape 604 pour maintenir une 
complexite du treillis d'estimation a un niveau, desire. 

L'elagage de l'etape 604 peut comporter diff6rentes realisations correspondant aux 
15 realisations de l'etape 504 de la figure 14. Les r6alisations peuvent etre combinees ou 
employ6es separement. 

Ainsi, dans une premiere realisation de l'elagage de l'etape 604, seuls les W noeuds les plus 
vraisemblables d'un 6tat (X», KJ sont gardes et ceci a partir d'un nombre n = x2 pour 
2 0 garantir la significativite des probabUites calcul6es (appelees vraisemblances). Ainsi, pour 
chaque etat (X„, KJ, seules les W probabilit6s les plus grandes parmi les probabilit6s 
calculees de decoder une sequence de K„ symboles connaissant des bits du flux de seconds 
elements binaires P(S 1 k /Yi Nk ) sont memorisees. 

25 Dans une variante, fl est possible de conserver seuls les w noeuds les plus vraisemblables 
pour chaque valeur possible de K, (correspondant a l'horloge bit). K„ est le nombre de 
symboles emis pour n bits recus lors du decodage. La relation entre W et w est donnee en 
(20) avec et K,,^ dependant respectivement de la probabiUte de la sequence S t k la plus 
et la moins vraisemblable. 
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Dans une deuxieme realisation de l'elagage de l'etape 604, le critere d'eiagage est base sur 
un seuil applique sur la vraisemblance des chemins. La vraisemblance d'un chemin est 
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donn6e par l'equation (13) (horloge bit). Dans cette equation, les termes P(Y|U) sont lies au 
modele de canal comme indique dans l'equation (16). Comme propose pour la figure 13, 
l'elagage consiste a supprimer tous les chemins pour lesquels la probabilite d'erreur binaire 
moyenne a la sortie du canal est superieure a p. Ainsi, a chaque instant bit n, le critere 
5 d'elagage comprend un seuil p compare aux probability P(Y n |U„) de maniere a ce que tous 
les chemins du treillis d'estimation pour lesquels est v6rifie P(Y n lU„) < (l-p) n , soient 
supprimes. Seules sont memorisees les probabilites de decoder un prochain symbole 
connaissant des bits du flux de seconds elements binaires P(S 1 Kn /Y 1 n ), calcul£es a partir des 
probabilites PC^JU,,) plus grandes que le seuil minimal p. Ainsi, a chaque instant bit n, tous 
1 0 les chemins du treillis d'estimation pour lesquels le critere d'61agage PC^JU,,) < (l-p) n est 
rempli, sont supprimes. 

Apres incrementation de la valeur n d'une unite, la probabilite P(S 1 Kn /Y 1 ^> pour le nombre 
de bits n (avec n>x2) est calculee a 1'etape 606. L'elagage realise a 1'etape 604 estreitere 
15 pour ce nouvel instant bit K^. Les 6tape 604 et 606 sont effectuees pour une valeur n 
increment6e jusqu'a la valeur N. Al'6tape 610, une fois la valeur N atteinte, la s6quence 
la plus vraisemblable est choisie, c'est-a-dire la sequence pour laquelle la probabilite 
P(S 1 K /Y 1 N ) est la plus grande parmi les sequences possibles du treillis d'estimation obtenu. 

20 A la fin de lestimation, les probabUitg P(S, K /Y 1 N ) marginales a posteriori sont connues pour 
toutes les sequences s x K restant apres 61agage. 

Les contraintes de terminaison mentionn6es peuvent 6tre vues comme un moyen de forcer 
la synchronisation a la fin de la sequence de symboles. En effet, elles contraignent le 
25 d6codeur a obtenir la bonne correspondance entre le nombre de. symboles decod6s et le 
nombre de bits observes. Cependant, ces contraintes n'ont pas d'effet de synchronisation sur 
le milieu de la sequence. 

L'addition a la sequence de symboles d'information complementaire permet de favoriser la 
3 0 re-synchronisation sur l'ensemble de la sequence. Ainsi, des symboles supplementaires, ou 
marqueurs, sont introduits a des positions connues dans la.sequence de symboles. Le choix 
de ces symboles est arbitraire, et se base sur le principe des techniques decrites dans (2, 4, 
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11, 5, 14). A titre d'exemple uniquement, un des symboles dans l'alphabet peut etre choisi 
ou le dernier symbole code pent etre repete. Ces marqueurs sont insures a position connue 
suivant l-horloge symbole. Par consequent, la position des bits supplementaires correspon- 
dants dans le train binaire d6pend de la sequence de symboles codee, et est aleatoire. 

Bien evidemment, les precedes developpes precedemment doivent tenir compte de cette 
information suppl6mentaire. Cette connaissance a priori peut etre exploitee par les procedes 
d'estimation. La variable K„(ka horloge symbole) indique quand un marqueur est attendu. 
Les probabilites de transition correspondantes dans le treillis ^estimation sont mises a jour 
en consequence. Une probabilite nulle est affectee a toutes les transitions qui n'emettent pas 
les marqueurs attendus, alors quWeprobabiUt6 de 1 est affect6e aux autres. Par consequent, 
certains chernins dans le treillis peuvent etre supprim6s, ce qui amene une reduction du 
nombre d'etats et une meilleure capacite de re-synchronisation. 

5 Des exemples de procedes de decodage robuste sont d6crits ci-dessus a titre d'exemple. 
D'autres variantes peuvent etre envisagees.De plus, d'autres criteresd'elagage peuvent etre 
pr6vues et combines avecles procedes de l'invention. La conversion M-aire/binaire n'estpas 
indispensable, mais permet d'utiliser de petites valeurs de T. Les procedes de decodage 
P r6sent6s ci-dessus pour des sources binaires peuvent 6galement fonctionner pour des 
0 sources M-aires. 

Comme developp6 dans la description ci-dessus, l'invention propose un procede de 
decodage robuste de codes arithm6tiques base sur l"utiUsation de codes de type codes 
arithmetiques comprenant les codes arithm^tiques a precision reduite, appel6s codes qua- 
2 5 si-arithmetiques. La precision de ces codes quasi-arithmetiques est un parametre ajustable 
en fonction de l'effet recherch6. Une precision tres elevee permet d'atteindre des performan- 
ces en compression identiques aux codes arithmetiques. Une precision plus faible permet 
d'accroitre la robustesse, au prix.dhme perte raisonnable en compression. 



3 0 Comparee a un decodage dit "instantane" classique, l'invention permet de reduire fortement 
leph6nomene de desynchronisation du d6codeur. Compar6e a d' autres m6thodes classiques 
de decodage dite "sequentielles" ou sur structures d'arbres, 1' invention offre une plus grande 
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souplesse pour ajuster un cx>mpromis entre Tefficacite de compression/la robustesse/la 
complexity du decodage. En effet, l'invention comprend un ensemble de precedes pouvant 
etre choisis selon la complexite peimise par le terminal, selon le taux d'erreur observ6 sur 
le canal. Un procede de decodage sequentiel comprend un mScanisme de synchronisation 
5 couple a une technique d'61agage qui permet de controler la complexite du processus de 
d6codage. Un proc6d6 de decodage sur treillis avec parametrisation du compromis 
complexite/robustesse/compression permet d'utiliser le procede de maniere couplee avecun 
decodeur de canal (decodeur de codes convolutifs), utilisation difficile avec les proc6des 
sequentiels traditionnels. 

.0 

Dans I'invention proposee, il doit Stre clair que les elements d'ordre mathematique decrits 
n'ont pas d'interet en eux-memes, mais seulement dans la mesure oil ils concernent des flux 
d'informations ( par exemple signal sonore, video), soumis a des contraintes materiellesvV 



1 Annexe 1 
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P{C k+1 \C k ) = P{S k+1 \C k ) (1) 



IJ P(S*|C fc _i) (2) 

fe=K„+l 



X = arg max P(X = as| Y) . (3) 

X 



N 

F(X,K\Y) = IJ nX n , K n \Y). (4) 
n=l 



Jr(.Jtn,-K-nK j — ' (5) 

OC P(X n , ifnllf ) - PO^i |X„, X„), 



P(X n = » n , i£ n = fc n |*i n ) = 13 P(X n _i = Sn-1, ifn-1 = fcn-l|lT -1 )- 

P(X„ = ap„, iC n = k n \X n -T. = X n -u Kn-l = K-lY 
P(Dn = W(x n _i,Jfe n _i)(»» > fc»)l 1 «)- 

P(X n+ i — Xn+x, K n+1 = k n+ i\X n = X n , K n — kn)- 
P(Cn+l = «(ah,,fe»)(*»+iJWl)l y »+ 1 )- 



p(5f|y^) = P(5f).p(t7fiy^) , 8 * 

oc P(Sf) -PCY^IC/f), 



P(5T*|Yf fc ) oc P^llf^J.PCfliklft-i). rq . 
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p«iL 1+1 i^_ 1+ i) 

P(Si = siir*^ 1 =y Nl=ni ) oc P(Si = «i|Ci)- 

II P C^ = yn|J7n = «n) if N k > N k -r 
1 sinon. 

W(S? n \Y?) oc P(Sf ») • p(y x n |i^) 



(10) 

(11) 

(12) 
(13) 



(14) 



p(^:_ 1+1 i5f- 1 ) = 

n p^iCfc-o si > jv- n _! 

1 sinon. 



(15) 



p(sf«|i7) oc p^f'-^-pfi^- 1 !^- 1 ). 



oc p(5f- i |y 1 ft ~ 1 )-ip(^«l«'»)- II i^sifcic5fc-i) 



(16) 



P(X n = x n ,K n = kn\Xn-l = X n - U K n -i = K^) = JJ F(S k \C k ^) (17) 

HYn\U n ) <(l~p) n (18) 

W = (JV^oa. - Nknin + l)u> (19) 

W = (Knmax ~ K nmin + l)w (20) 
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P(Yi = Q\Ui = 0) = 1 
P(Yi = 1|CT« = 0) 
¥(Xi = OlDi = 1) 

HYi = Wi = i) = i 
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Rcvendications 

1. Proced6 de decodage comprenant les etapes suivantes: 

a. prevoir un modele source 

5 propre & definir des probability pour des symboles x-aires (S^, 

associes a des transitions entre etats sources (C*) dans un diagramme d'etats source, 
un modele dit "produif 

propre & definir une correspondance entre les symboles x-aires et de premiers Elements 
binaires, 

10 en fonction d'etats produits dans un diagramme d'etats "produit", un etat produit 

etant fonction d'un etat source, 

b. recevoir des propriety de "canal", 

propres h definir des probability de recevoir de seconds elements binaires, sur emission 
de premiers elements binaires, 
15 c. pour un flux code re<ju de seconds Elements binaires, 

calculer des probability de decoder des symboles x-aires, 
connaissant des bits du flux de seconds elements binaires, 

ou des probability d'aboutir dans des etats produits connaissant des bits du flux de 

seconds elements binaires, 
20 a partir des proprietes de canal, du modele source et du module produit, 

d. reconstituer le flux de symboles x-aires le plus probable a partir des probability calculees 
i Petape c. 

2. Precede selon la revendication 1, caracterisee en ce que le modele produit est etabli a 
25 partir du modele source et d'un modele de transfert, propre a etablir une correspondance 

entre des symboles x-aires et de premiers elements binaires, en fonction de probability 
assoctees k des symboles x-aires. 
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3. ProcSde selon la revendication 2, caract6rise en ce que le mod&le source comprend des 
transitions associfies & une correspondance entre des symboles m-aires (A*) sources et des 
symboles x-aires cibles (SJ, en particulier m^x et x=2. 
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4. Precede selon l'une des revendications precedentes, caracterisee en ce que les modeles 
de produit et de transfert sont de type arithmetique et en ce que le modele de transfert est un 
modele de codage ou de decodage. 

5 

5. Proced6 selon Tune des revendications precedentes, caracterise en ce que l'6tape c. 
comprend le calcul, par etapes, des probabilites de d6coder un symbole x-aire connaissant 
des bits du flux de seconds 616ments binaires. 

10 6. Precede selon l'une des revendications pr6cedentes, caracterise en ce que l'etape c. 
comprend, pour chaque etape a partir d'un nombre donn6 d'etapes et selon un critere de 
selection donn6, la selection de certaines probabilites parmi les probabilit6s de decoder un 
symbole x-aire, connaissant des bits du flux de seconds elements binaires, calculees lors 
d'une etape. 

15 

7. Precede selon la revindication 6, caract6rise en ce que le critere de selection de l'etape 
c. comprend un nombre fix6 de probabilites a garder parmi les plus grandes des probabilites. 

8. Precede selon l'une des revendications 6 et 7, caracterise en ce que le critere de selection 
20 de l'etape c. comprend un seuil minimal compar6 aux probabilites de recevoir de seconds 

elements binaires, sur emission de premiers elements binaires, de maniere a ne memoriser, 
a chaque etape, que les probabilites de d6coder un symbole x-aire, connaissant des bits du 
flux de seconds elements binaires, calcul6es a partir des probabilites de recevoir de seconds 
elements binaires, sur emission de premiers elements binaires, plus grandes que ce seuil 

2 5 minimal. 

9. Precede selon l'une des revendications 5 a 8, caracterise en ce que l'6tape c. comprend 
l'elaboration d'un arbre constitue d'etats reliees par des transitions, chaque etat correspon- 
dant a une etape de calcul de probabilites et chaque transition, au depart d'un etat, 

3 0 correspondant a 1'une de ces probabilit6s calculees. 
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10. Procede selon les revendications 1 et 9, caracterise en ce que Fetape d. comprend le 
calcul du produit des probabilites correspondant a des transitions successives dans Farbre 
et la selection du produit le plus grand correspondant aux transitions successives les plus 
probables et au flux de symboles x-aires le plus probable. 

11. Procede selon les revendications 5 et 10, caracterise en ce que P6tape c. comprend le 
calcul, pour un nombre donne de symboles x-aires successif s a decoder, des probabilites de 
decoder une sequence de symboles x-aire connaissant des bits du flux de seconds elements 
binaires. 

12. Procede selon les revendications 5 et 10, caracteris6 en ce que Fetape c. comprend le 
calcul, pour un nombre donne de bits successifs du flux de seconds elements binaires, des 
probabilites de decoder une sequence de symboles x-aire, connaissant des bits du flux de 
seconds Elements binaires, 

13. Procede selon Tune des revendications 1 h 4, caracterise en ce qu'& partir du modele 
produit etablissant ladite correspondance, Petape c. comprend 

cl. la construction successive d'etats relies par des transitions d'un diagramme d'etats 
produit depuis un etat d'origine a un etat final, chaque etat definissant un nombre attendu 
de symboles x-aires decodes pour un nombre donne de bits du flux de seconds el6ments 
binaires depuis P6tat d'origine. 

14. Procede selon la revendication 13, caracterise en ce que la construction de Tetape c. 
comprend 

c2. la suppression successive d'etats du diagramme d'etats produit, depuis Petat final a Tetat 
d'origine, pour lesquels le nombre attendu de symboles x-aires pour un nombre donne de 
bits du flux de seconds elements binaires est different d'un nombre connu maximal de 
symboles x-aires. 

15. Proc6d6 selon Yxme des revendications 13 a 14, caracterise en ce que Petape c. comprend 
cl. le calcul, apres chaque etat construit, des premieres probabilites successives d'aboutir 
dans un 6tat donne connaissant des bits passes du flux de seconds elements binaires, 
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c2. le calcul, pour chaque etat depuis Petat final a Petat d'origine, des secondes probability 
successives d'observer des prochains bits du flux de seconds elements binaires connaissant 
cet etat donne. 

5 16. Procede selon Tune des revendications 13 a 14, caracteris6 en ce que Petape c. comprend 
c3. le calcul, pour chaque etat du diagramme depuis Petat d'origine k Petat final, des 
premieres probabilites successives d'aboutir dans cet etat donne connaissant les bits passes 
du flux de seconds elements binaires, et , pour chaque etat depuis Petat final a Petat 
d'origine, des secondes probabilites successives d'observer des prochains bits du flux de 
1 0 seconds elements binaires connaissant un 6tat donne. 

17. Procede selon Pune des revendications 13 a 16, caracterise en ce que Petape c. comprend 
le calcul de la probability d'aboutir dans cet etat donne connaissant les bits du flux de 
seconds 61ements binaires en effectuant le produit des premieres et deuxiemes probabilites 

1 5 pour chaque etat. 

18. Procede selon les revendications 1 et 17, caracterise en ce que Petape d. comprend le 
calcul des produits de probabilites des differents etats successifs possibles puis la s61ection 
du produit maximal. 

20 

19. Procede selon les revendications 13 et 18, caracterise en ce que Petape d. comprend 
Petablissement du flux de symboles x-aires correspondant a ces etats successifs en utilisant 
la correspondance du modele produit. 

2 5 20. Decodeur comprenant 

une premiere entree pour recevoir un flux cod6 de seconds elements binaires, 
caracterise en ce qu'il comprend 

une deuxieme entree pour recevoir des proprietes de canal propres a definir des probabilites 
de recevoir de seconds elements binaires, sur emission de premiers elements binaires, 
30 un premier module definissant 
- un modele source 
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propre a definir des probabilites pour des symboles x-aires associes & des transitions entre 
etats sources dans un diagramme d' etats source, 
- un modele dit "produit", 

propre a d6finir une correspondance entre les symboles x-aires et de premiers Elements 
5 binaires, 

en fonction d'etats produits dans un diagramme d' etats "produif % un etat produit 
6tant fonction d'un etat source, 
un module de traitement 

adapt6 pour calculer des probabilit6s de d6coder des symboles x-aires, 
1 0 connaissant des bits du flux code de seconds elements binaires, a partir des 

proprietes du canal, du modele source et du modele produit, 
ou des probabilit6s d'aboutir dans des etats produits connaissant des bits du flux de 
seconds elements binaires, 
et adapte pour reconstituer le flux de symboles x-aires le plus probable si partir de des 
15 probabilites calculees. 
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